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Polariton waveguides from a quantum dot chain in a nanowire photonic crystal: 
an architecture for quantum information science and waveguide QED 
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We introduce a polariton-waveguide structure, comprised of a nanowire-based photonic crystal 
waveguide with a quantum dot embedded in each unit cell. Using realistic designs and parameters, 
we derive and calculate the fundamental electromagnetic properties of these polariton waveguides, 
with an emphasis on the photon Green function and local optical density of states (LDOS). Both 
infinite and finite-size waveguides are considered, where the latter’s properties are calculated using 
a Dyson equation approach without any approximations. We demonstrate dramatic increases, and 
rich fundamental control, of the LDOS due to strong light-matter interactions in each unit cell 
through periodic quantum dot interactions. Consequently, these structures allow the exploration of 
new regimes of waveguide quantum electrodynamics. As an example application, we consider the 
coupling of an external target quantum dot with a finite-sized polariton waveguide, and show that 
the single quantum dot strong coupling regime is easily accessible, even for modest dipole strengths. 

PACS numbers: 42.50.Ct, 42.50.Nn, 78.67.He, 78.67.Qa 


The design of nanoscale solid state devices to control 
light-matter interactions is highly desirable, both for ap¬ 
plications in quantum information science [1] and to ex¬ 
plore novel regimes of quantum electrodynamics (QED) 
[HIS]. Photonic crystals (PCs) with embedded quantum 
dots (QDs) have had much success in this regard by mod¬ 
ifying the LDOS [1H3, although they have been hindered 
by fabrication issues including surface roughness [MIQ] 
and limited control of the QD properties HU. Arrays of 
nanowires (NWs) grown through a molecular beam epi¬ 
taxy (MBE) technique HU [T3] have been proposed as 
an alternative PC platform to help mitigate these issues 
HU HU , with the potential to contain identical or very 
similar QDs inside each NW of a given radius HU HU- 
This opens the idea of NW PC waveguides where each 
waveguide channel NW contains an identical QD embed¬ 
ded in its center. Coupled dipole chains in free space have 
interesting collective properties, and can act as subwave¬ 
length waveguides |18j : implementing such chains in PC 
waveguides could result in new physical behavior or im¬ 
proved performance. Metamaterials, comprised of multi¬ 
ple elements that are engineered to have unique and useful 
properties, have been used to produce exotic waveguides 
with a dramatically different LDOS than simple dielectric 
structures, e.g., manifesting in large spontaneous emission 
rate enhancements HU US] for dipole emitters, even with 
material losses. Plasmonic polariton waveguide structures 
m and metal nanoparticle chains coupled to traditional 
waveguides [22] have demonstrated chain-mediated cou¬ 
pling between plasmonic and photonic excitations leading 
to an anti-crossing in the band structure. This polari¬ 
ton anti-crossing is also observed with quantum wells in 
ID distributed Bragg reflector waveguides [53]. All these 
structures essentially lead to normal mode splitting. 

Inspired by the metamaterial concept and advances in 
NW growth techniques, the inclusion of a QD chain in a 
PC waveguide can similarly reshape and greatly enhance 


the system LDOS and result in drastically different prop¬ 
erties and behavior, but without the large metallic losses 
of plasmonic metamaterials. In this Letter, we introduce 
and explore the optical properties of a nanophotonics sys¬ 
tem comprised of a PC waveguide with an embedded 
QD in each unit cell, which we describe as a “polari¬ 
ton waveguide” because its excitations are mixed light- 
matter states due to the collective coupling of the QDs 
with the waveguide Bloch mode. Eigure [^a) shows a 
schematic of our proposed structure. To elucidate the 
underlying physics of these systems, we derive the pho¬ 
ton Green function (GF) of both infinite and finite-sized 
systems polariton waveguides, and explore their coupling 
with a single external “target” QD, which is found to be 



FIG. 1. (Golor Online) (a) Proposed polariton waveguide, 
with yellow QDs embedded inside the PG waveguide channel, 
and a red external QD coupling to the structure. Red arrows 
denote the waveguide direction and NWs are cut out so em¬ 
bedded QDs can be seen, (b) Slow light region of the band 
structure on top, and Im{G(r„,r„)} in units of (see text) 
on bottom; the infinite polariton waveguide in solid blue is 
compared with the original PC waveguide in dashed orange. 
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in the strong coupling regime. This achievement of strong 
coupling between a single isolated emitter and a waveg¬ 
uide mode, to the best of our knowledge first reported 
here, allows for the creation of devices relying on quan¬ 
tum cavity physics, such as photon blockades and single 
photon switches 0 0, with reliable input/output cou¬ 
pling on chip. Our polariton waveguide architecture can 
also be adapted for other systems such as circuit QED [2] . 

Our proposed structure exploits the elevated NW PC 
waveguide design of Ref. [T3J where a PC waveguide is 
formed from an organized array of GaAs NWs extended 
from an AlO substrate [23], with the waveguide chan¬ 
nel introduced by reducing the radius of a row of NWs 
from rt = 0.180a to = 0.140a along e^,. Importantly, 
this structure is based on current fabrication techniques 
and properties are determined through full 3D calcula¬ 
tions including radiative coupling effects. The PC waveg¬ 
uide with lattice constant a = 0.5526 pm contains a single 
vertically-polarized below-light-line waveguide band with 
a mode edge near the telecom wavelength of 1.550/rm 
(ranging from 755 —795meV). Light is confined to the up¬ 
per GaAs portion of the NWs (height 2.27 a), while the 2 a 
AlO layer separates them from the substrate and an array 
width of 7 a is sufficient to prevent in-plane losses. As this 
waveguide band approaches the mode edge it flattens due 
to symmetry and the group velocity goes to zero, causing 


the LDOS to diverge as seen in Fig. 1 (b) an effect that is 


inevitably spoiled by disorder-induced losses in real sys¬ 
tems |8]. To form the polariton waveguide of Fig. i[a) 


we embed a QD in the center of the GaAs layer of each 
waveguide NW at r„ g = rg-l-naej;, where n is an integer. 
We take the embedded QDs to have a Lorenzian polariz¬ 
ability a = a{oj)ez = 2 wo|dpe 2 /(Iieo(wg — — iPow)), 

with a dipole moment |d| = SOD (0.626e-nm) and po¬ 
larization decay rate Pg = l/reV (including both non- 
radiative decay and coupling into non-waveguide modes); 
these are similar to experimental parameters for InAs 
QDs at 4K [2S|. The QD exciton line is taken to be 
ujQ = 794.5 meV, which is in the moderately-slow-light 
regime of the waveguide band resulting in a group veloc¬ 
ity Vg = c/30.4 and = 40.1 (relative LDOS enhance¬ 
ment, Fz = Im{G}/Im{G*'}). Throughout this work, we 
connect to the GF G(r, r/w), which describes the sys¬ 
tem electric-field response at r to a point source at r' 
(fully including all light scattering events), and whose 
imaginary and real components at equal space points are 
directly proportional to the system LDOS (and emitter 
spontaneous emission rate), and Lamb shift, respectively 
[DUS]. These GFs are projected along the relevant 
component for vertically polarized QDs: G = e^, • G • 
and given in units of the imaginary part of free-space GF 
Im{G*'(r, r;a;)} = w^/(67rc^) = p^{uj) [2S]- When we con¬ 
sider the addition of an external target QD to the waveg¬ 
uide, this will be embedded on the top of waveguide NW 
where n indexes the unit cell, as seen in 


Fig. 1 'a) both for fabrication purposes and because this 
is where the antinode of the waveguide mode resides. 


We first consider an infinite embedded QD array to 
help explain these structures’ underlying physics, exploit¬ 
ing Bloch’s theorem and treating the QDs as a pertur¬ 
bation to each PG unit cell Ae(r) = (5(r — rg)Q;(a;(,)e 2 , 
which shifts the waveguide resonance cok for a given 
k = kx to We obtain from perturbation theory, 

Ae(r)|e 2 •Ufc(r)pdr) [27H21], where 14 is 
the unit cell volume and Ufc(r) is the waveguide unit-cell 
function, normalized through fy e(r)|ufc(r)pdr = 1, with 
corresponding waveguide mode ffe(r) = y^Ufc(r)e*^’^ 
where L the waveguide length. The waveguide band is 
split by the QD-Bloch-mode interaction, resulting in a 
pair of complex eigenfrequencies w/ = uJk,± + irk,±/2 


at each k, with ujk,± = ± f 


wg - cjk) + 4g^ 


and 


- k,± 


= Tn 




The QD-unit-cell cou- 


pling parameter g/- = • Ufe(rg) gives the strength 

of the light-matter interaction leading to anti-crossing and 
is the continuum form of g, the single mode quantum op¬ 
tical coupling rate. This modified band structure is com¬ 
pared with that of the original PG in Fig. 1 where 
the strong mode splitting (g^g = 120 peV Pg ) flat¬ 
tens the dispersion as (jJk,± ojo and losses remain low: 
L/c.i < Lg ^ 0Jk,±. The inclusion of the QD array has 
caused the waveguide band to split, corresponding to a 
mixed light-matter excitation: the polariton waveguide. 

The above approach was also applied to the waveguide 
Bloch modes: away from uJk,± —> wg where perturbation 
theory naturally breaks down we hnd only the waveguide 
Bloch mode contribution is signihcant and Ufc ± = u^, 
and we only show results for frequencies where perturba¬ 
tion theory still holds. We write the GF of this inhnite 
polariton waveguide as a sum over these waveguide Bloch 


modes, Gp(r, r'; w) = 




and con- 


LO'±{k)'^—ljj'^ 

verting this sum to an integral in the complex plane [29] , 
we arrive at an analytic expression similar in form to that 
for a regular PG waveguide [11150]: 

Gp(r, r'; cu) loix - x')uk^ (r)u*,^ 

2Vg L 


-b0(a:' - a;)ut (r)ufc_^(r')e 




, ( 1 ) 


but with Vg(uj) = Vg(uj) - ^ I k^, where Vg(uj) > ^ 

is the polariton waveguide group velocity and Kg, = 
emphasize that both Vg(uj) and k^, are 
found from the polariton waveguide oj± — k relationship, 
which differs greatly from that of the original PG near wg 
due to the polariton splitting. This was derived assum¬ 
ing Kgj kgj, which remains valid until v'g ~ c/2300, 
where |w — cugl < 20peV. The effect of the QD ar¬ 
ray, demonstrated in Fig. is thus to produce large 
and tunable LDOS enhancements near resonance by flat¬ 
tening the band structure; as one moves away from wg, 
the original PG waveguide GF is quickly recovered. At 
the above minimum detuning, these polariton waveguides 



































3 



FIG. 2. (Color online) (a) Real and imaginary compo¬ 
nents of , r 5 i;u;) in solid blue, compared with 

G(®i)(r26 , r 2 e',ou) in dashed orange and = 0 results in dash- 
dotted red. (b) |G'^“^^(r„ , r 5 i;aj)| for N = 101 polariton 
waveguide near ojq, where An = n — 51 and the array ends 
at An = ±50. All values are in units of 


yield = 3100—a dramatic rate enhancement and in¬ 
dicative that these structures should enter the strong cou¬ 
pling regime of QED, even at the single quantum level. 

We next consider polariton waveguides with a finite 
number of embedded QDs, more representative of real 
systems. We can no longer exploit Bloch’s theorem, but 
instead include QDs iteratively in the system GF via an 
exact Dyson equation approach. Denoting the background 
PC waveguide Green function as and starting with 
unit cell n = 1, we introduce a QD at r„ 0 ; calculate 
the resultant Green function including scattering 

from this QD, and then use this as the background Green 
function to introduce a QD in the subsequent unit cell. 
From the Dyson equation, QD n can be included self- 
consistently in the system GF through m G(")(r,r') = 
G("-i) (r, r') ± G(”-i) (r, r„,o) • a • G^”) (r„,o, r')• This was 
done numerically using the same system as for the infinite 
case (i.e. identical PC waveguide and QDs) for QD chains 
of increasing length, TV = 1 —>■ 101. 

Figure shows corresponding to a chain suffi¬ 

ciently long to produce substantial LDOS enhancements 
and begin to recover the infinite chain result, while still 
demonstrating important finite-size effects. We chose to 
emphasize the position on top of the central NW of the 
QD array rsi (where a target QD will be later embed¬ 
ded) because constructive interference from repeated QD 
scattering maximizes the polariton effects. The addition 
of a single QD introduces a dip in G at wq, however with 
increasing N the build up of off-resonant enhancement 
from QD scattering leads to the formation of a strong 
resonance in Im{G} which exceeds G*^°^ for TV > 15. This 
resonance is red-shifted from ojq, with a peak = 338.14 
and FWHM Pi = 3.35/reV at wi = luq — 21.49/reV 
for the 101 QD case. As more QDs are added, this 
polariton peak grows, narrows, and blue-shifts towards 
Wo while additional weaker red-shifted resonances begin 
to appear. Above wq, resonances also form which grow 


and become more numerous with increasing N, indicating 
that they arise from QD chain Fabry-Perot (FP) modes. 
The higher two FP modes for the N = 101 structure are 
at wpp/ = wi ± 25.12/reV and wpp = Wi ± 29.02 ^eV. 
Away from wq, converges to the PG waveguide GF, 
as was seen for the infinite case. The large number of res¬ 
onances in these waveguides also produce a richly varying 
Re{G} which remains substantial over a broad frequency 
region, particularly for larger N. 

The above results were also found to be robust to 
disorder in embedded QD position and dipole moment; 
e.g., random variations of up to 10nm and 3D reduced 
the strength of the primary peak slightly (~ 2%) and 
introduced a frequency shift of ~ 1 p.eY, but the un¬ 
derlying system behavior remained unchanged. Simi¬ 
larly, these findings were also verihed to hold for Tg = 
0—10 /ieV, with the wg dip extending deeper and res¬ 
onances initially appearing at lower N and closer to 
Wg as Tg is reduced. Figure |^Kb)| shows the propa¬ 
gator |G(i°i)(r„ , r 5 i;w)|, which is proportional to the 
coupling strength between the target QD location and 
various points in the structure, and far greater cou¬ 
pling rates are found than for a bare ideal PG waveg¬ 
uide. For instance, |Im{G^^°^^(r„,r 5 i;w)}| > 300p^(w) 
and |Re{G(^°^)(r„,rgi;w)}| > 150p^(w) are found for 
An < 15, and one can produce effectively any arbitrary 
combination of Im{G} and Re{G} through careful choice 
of separation and frequency. This GF reshaping persists 
even past the mode edge of the polariton waveguide: once 
one is outside of the QD chain, G(^)(r„,r„/;w) experi¬ 
ences only a phase shift of I as n and n' are 

varied and |G^^°^^(r„,r 5 i;w)| has a peak of 101.8p^(w) 
at « wi (wg — 21.6/reV), more than double the bare PG 
waveguide result: |G(°^(r„,r„/;wi)| = 49.5p^(wi). No¬ 
tably these finite-sized GFs are derived without any ap¬ 
proximations, and reproduce the physics of the infinite 
structure, namely strong light-matter interactions result¬ 
ing in dramatic enhancements in the system LDOS. 

We now consider an important quantum optical ap¬ 
plication of these polariton waveguides, studying the 
interaction of a single external QD at r* = rgi of 
the N = 101 polariton waveguide. The LDOS en¬ 
hancements, specihcally the polariton peak at wi, are 
possibly sufficient to strongly couple the waveguide to 
a realistic QD. We follow a first principals quantiza¬ 
tion procedure appropriate for lossy inhomogeneous sys¬ 
tems [35] such as these polariton waveguides. The sys¬ 
tem Hamiltonian in the dipole approximation, consist¬ 
ing of a single target QD interacting with the elec¬ 
tromagnetic environment of the polariton waveguide, is 
given hy H = /iwt(j+(j“ + f dr dcJi ■ 

b(r;wi) — ((T+ ± i7“)(dt • E(rt) + H.c.), where the QD 
is a two-level atom (TLA) with exciton frequency Wt 
and transition dipole moment d^; b(r;w;) is the bosonic 
field annihilation operator associated with the field fre¬ 
quency uji which generates the electric-held operator via 
E(r)oc/ duJi f dr'^Im{e(r';uJi)}G(r, r'; w;) •b(V; w/) [32]. 
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(a) (b) 



FIG. 3. (Color online) S^{uj) in dashed blue and S{vd,u!) in 
solid orange for ujt = (jJi, with dt = 30 D in (a) and cdt = Wfp 
with dt = 60D in (b). Ini{G(rf, rt; a;)} in arbitrary units is 
shown in dash-dotted gray. 



We Laplace transform the resulting operator Heisen¬ 
berg equations of motion to construct the spontaneous 
emission spectrum at detector position ru: S{ru,uj) = 
(El(r£i,a;)E(rD,a;)), which yields 


S{rD,u}) 


{ujt -ba;)G(rD,rt;a;) ■ d/ep ^ 
lJI — up- — a;E(a;) — iioTt 


( 2 ) 


where we have assumed an initially excited QD in vac¬ 
uum. The self-energy is E(w) = ^ 

is the polarization decay rate of the target TLA due to 
interactions with the environment (i.e., anything other 
than the LDOS given by the polariton waveguide GF). 
The spectrum depends not only on the GF at TLA 
position r* but also on the propagator to the detector 
G(r£),r*; w). When we present spectra below, we show 
both S{r£),uj) and where S°{uj) is calculated from 

Eq. ([^ by replacing the propagator in the numerator with 
This is done because S{rD,uj) is often dis¬ 
torted by the rich spectral features of the GFs of these 
waveguide structures (see Fig. 2b). The bare spectrum 
S^{u!) thus more cleanly shows the system energy levels 
and dynamics of the target QD. 

We assume a e^-aligned target QD with a polariza¬ 
tion decay rate F^ = 1/ieV and calculate spectra for tar¬ 
get QD dipole moments of 10, 30, and 60 D, encompass¬ 
ing the range of QDs which could feasibly be coupled to 
the polariton waveguide. For a simple Lorentzian LDOS, 
the resultant system dynamics can be approximated with 
the Jaynes-Cummings (JC) Hamiltonian, where the in¬ 
teraction of a single photonic mode and TLA leads to 
an anti-crossing as they approach resonance (splitting of 
2hg on resonance); near wi we define an effective coupling 


constant ^ 

target dipole moments produce ges of 1.24, 3.72, and 
7.43 /reV, respectively, which meet the criteria for strongly 
coupling with the Wi resonance: 2geg > Fi,Fi [33]. We 
assign a detector position td = r* — 55aea;, which is 
outside the polariton waveguide portion of the structure 
to reduce the filtering of the detected spectrum, while in 
the same position in the unit cell as the QD to maximize 
coupling. The emitted and detected spectra S'°(a;) and 


FIG. 4. (Golor online) Peaks in S^{u!) in blue as a function of 
ojt and dt- Results for dt = 30 and 60D are presented on the 
bottom and top, with uit near cui and in the FP region on the 
left and right, respectively. Crosses and circles denote peaks 
of Figs. and |(b)[ wt is in dashed red, and LDOS peaks 
are in dashed gray. 


S{rD,u}) at Wt = uJi for the 30 D QD is shown in Fig.|i]Ka) 


Remarkably, substantial splitting is seen in both oPuj) 
and 5'(rD,a;) despite the modest choice of dt and peak 
locations of a;_ = uii — 3.02 /reV and ui+ = uji+ 3.67 /xeV 
agree well with the JC Hamiltonian when the Lamb shift 
is included; however we stress that the widths and weight¬ 
ing of the spectral peaks, as well as propagation effects 
and polarization decay, are only captured through the 
formalism leading to Eq. in particular, the reduced 
height of the uj- peak is a result of the unique shape of 
Re{G} near wi producing a large positive Lamb shift. 
The importance of propagation effects can clearly be seen 
by comparing 5'°(a;) and iS'(r£),w), where the depletion 
above oJi reduces the height of the a;+ peak. In Fig. ^’b) 


we show the spontaneous emission spectra at Wj = wpp 
for dt = 60 D, where the QD-polariton-waveguide cou¬ 
pling is sufficiently strong that significant exchange oc¬ 
curs with both the wpp and wpp' modes for wt = wpp. 
This results in a triplet forming in S^{uj), with peaks at 
wpp' — 0.44^eV, wpp — 1.81/ieV, and wpp -|- 2.16^eV; 
the splitting is substantially stronger than rFp/2 so these 
peaks are clearly seen at the detector position as well. 


Finally, to demonstrate that the features in the above 
spontaneous emission spectra are a consequence of the 
strong coupling regime, the clearly resolvable peaks of 
S^{u!) as uJt is brought to resonance with wi are shown 
in Fig. [^ for a target QD with dt = 30 and 60 D. The 
left plots depict anti-crossing with the primary resonance 
at wi, while the peaks as Ut is swept through the FP 
region are shown on the right and markers denote the 
peaks of Fig. [^ Remarkably for a waveguide system, 
all three QDs {dt = lOD is not shown) demonstrate a 
clear anti-crossing as they approach wi, conclusive evi¬ 
dence that they are strongly coupled to the polariton- 
waveguide despite the complicated nature of this system. 
While similar strong coupling behavior has been theoret- 
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ically predicted for coupled-cavity PC waveguides [33], a 
later study showed this was inevitably spoiled by disor¬ 
der [To] . Strong coupling has recently been observed with 
Anderson-localized cavities in disordered PC waveguides 
[5], although this is fundamentally different from strong 
coupling with a propagating waveguide mode as reported 
here. It is also impressive that our predicted splitting is 
seen even for the 10 D target QD, which is much weaker 
than that used in other studies. Of further interest is the 
system behavior in the FP region of the polarition waveg¬ 
uide, where the unusual shape of the system LDOS leads 
to multiple anti-crossings which are poorly described by 
a JC Hamiltonian. The interaction of the FP modes with 
the 60 D target QD is particularly striking: up to four 
energy levels are seen for a given ujt and the multi-mode 
nature of the QD-polariton waveguide has flattened these 
anti-crossing lines. 

In conclusion, we have proposed a new nano-engineered 
metamaterial system, a polaritonic waveguide, consisting 
of a PC waveguide with periodic embedded QDs. We de¬ 
veloped two separate approaches to describe the physics 
of this system, studying both infinite and finite polari- 
ton waveguides and demonstrating that in both instances 
strong light-matter interactions lead to rich and dramatic 
LDOS enhancements which are robust to disorder and 
without the losses typically associated with metallic meta¬ 
materials. We then considered the interaction of a finite- 
size structure with an external QD and we showed that 
these LDOS enhancements can be exploited to strongly 
couple with a single external emitter and produce inter¬ 
esting spectral features clearly visible in the externally 
detected spectrum. These structures could thus be used 
to design complex devices for quantum information and 
explore new regimes of open system waveguide QED. 

This work was supported by the Natural Sciences and 
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University. 
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